Cyclotron resonance in bilayer graphene 
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We present the first measurements of cyclotron resonance of electrons and holes in bilayer graphene. In 
magnetic fields up to B = 18 T we observe four distinct intraband transitions in both the conduction and 
valence bands. The transition energies are roughly linear in B between the lowest Landau levels, whereas they 
follow a/B for the higher transitions. This highly unusual behavior represents a change from a parabolic to a 
linear energy dispersion. The density of states derived from our data generally agrees with the existing lowest 
order tight binding calculation for bilayer graphene. However in comparing data to theory, a single set of fitting 
parameters fails to describe the experimental results. 

PACS numbers: 78.66.Tr; 71.70.Di; 76.40.+b.; 78.30.Na 
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Experiments into the properties of single layer graphene 
have demonstrated the existence of unusual charge carriers, 
analogous to massless, chiral Dirac particles having a Berry's 
phase of it and resulting in a new, half-integer quantum Hall 
effect |H 0]. Bilayer graphene has shown equally exciting 
properties, creating a system of massive chiral particles with 
a Berry's phase of 2ir and exhibiting another, distinct quan- 
tum Hall effect H. Both phenomena derive from the unusual 
energy dispersion of graphene leading to the presence of a dis- 



the presence 



tinctive Landau level (LL) at zero energy | 

The peculiar dispersion relations of graphene create unique 
LL spectra in the presence of a magnetic field, B. In single 
layer graphene the linear dispersion leads to distinctive y/B- 
and index-dependent energies for the LLs J3, Htl which have 
recently been observed in infrared (IR) cyclotron resonance 
(CR) ISlllOll- In contrast, bilayer graphene exhibits a very un- 
usual hyperbolic dispersion, which is predicted to have a more 
intricate index dependence and, uniquely, a varying B-field 
dependence. In a nearest-neighbor, tight-binding approxima- 
tion at zero magnetic field, bilayer graphene is a semimetal 
with degeneracies occurring at two inequivalent points at the 
corners of the hexagonal Brillouin zone 0. 11]. The low- 
energy part of the conduction and valence bands consists of 
two concentric hyperbolae, offset by the interlayer coupling 
energy 71, with energies given by 
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where c is the band velocity and h is Planck's constant. The 
lower energy branches meet at the charge neutrality point as 
depicted in the upper right inset to Fig. 1. This zero field 
dispersion has been explored via angle resolved photoemis- 
sion spectrometry (ARPES), demonstrating such a band ar- 
rangement ifTitl . In a magnetic field, the lower branches are 
expected to quantize into LLs with energies lfl3l [3 
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Here, the index n = ±1, ±2... describes a set of electron (+) 
and hole (— ) LLs which are four-fold degenerate, accounting 
for spin and valley degrees of freedom. The unusual zero- 
energy n = LL is 8-fold degenerate J3[0]. The term A is 
short for V^eBc^h, where e is the electron charge. With 71 
set to zero, Eq. 2 contracts to give the LL energies for single 
layer graphene. 

The B-field dependence of Eq. 2 differs sharply from both 
the LL energies of single layer graphene, and from traditional 
2D systems with a parabolic dispersion for which the energies 
are evenly spaced and linear in B. Instead, at low energies, 
E n <C 71, Eq. 2 approaches a linear-in-B spectrum, while 
at higher energies, E n ~ 71, it approaches a "\/B-dependence 



I1J, ll5|]. Thus the LL spectrum of bilayers is expected to 



move from that characteristic of a parabolic dispersion to that 
of a linear dispersion. 

Here we report CR measurements via IR spectroscopy on 
bilayer graphene, in magnetic fields up to B = 18 T. In both 
the conduction and valence bands, we resolve four LL transi- 
tions whose resonance energies decrease with increasing LL 
index and decreasing B-field. The B dependence of the res- 
onance energies is close to \/B for all data, except the tran- 
sitions between the lowest LLs for which the dependence is 
roughly linear in B. The observed energies agree qualitatively 
with the predicted change in B-field dependence, but we find 
the changeover to a \/B behavior to occur at lower energies, 
and more suddenly, than expected. 

The CR data can be exploited to arrive directly at the den- 
sity of states (DOS) for bilayer graphene, finding it to be 
roughly linear in energy with a constant offset, as expected for 
a hyperbolic dispersion. Details of this data reduction indicate 
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contributions to the CR energies that go beyond a simple, one 
particle band structure picture. 

The device used in this study consists of a 400 /im 2 flake 
of bilayer graphene extracted from bulk Kish graphite and 
deposited on a Si/SiC>2 substrate. This substrate is semi- 
transparent in the IR yet sufficiently conductive to serve as 
a gate, and is thinned and wedged to suppress Fabry-Perot in- 
terference of the incident IR light. Standard microlithography 
techniques are used to electrically contact the bilayer. The 
sample resides on a moveable stage at the focus of a parabolic 
cone. Transport properties are measured with standard lock-in 
techniques. All experiments are performed in liquid helium in 
a magnetic field perpendicular to the graphene sheet, utilizing 
a Bruker IFS 66v/S Fourier-transform IR spectrometer cou- 
pled with a composite Si bolometer. A gate voltage V g applied 
to the substrate varies the carrier density in graphene continu- 
ously from holes to electrons. Varying V g amounts to shifting 
the Fermi level through the LLs, changing the LL filling factor 
v and creating quantum Hall plateaus with R xy = h/ {ye 2 ). A 
representative trace of R xy taken at B = 18 T is shown as the 
upper left inset to Fig. 1 . Clear quantum Hall plateaus appear 
at v = ±4, ±8, and —12, demonstrating the 8-fold degener- 
acy of the n = LL and the 4-fold degeneracy for all other 
LLs, identifying this sample as a graphene bilayer OLI^]. 

The IR absorption of bilayer graphene is determined as in 
Ref. Time-averaged IR transmission scans are taken as 
a function of energy at constant magnetic field for two differ- 
ent filling factors ( V g values), and one is divided by the other. 
As a result, a transmission minimum in the numerator (de- 
nominator) appears as a dip (peak) on a background of unity. 
In all cases, we use v = ±4, where the Fermi level lies be- 
tween the n = and n = ±1 LLs, as the normalizing scan 
for all other transitions. Figure 1 shows three representative 
normalized traces at B -fields of 18 T, 14 T, and 10 T, with V g 
set so the corresponding filling factors are v = —8, +12, and 
+16 respectively. Most sharp noise spikes are due to harmon- 
ics of 60 Hz, while the noisy band around 90 meV is due to 
strong IR absorption in the optical window. For all traces, two 
broad resonances consisting of one peak (due to the normaliz- 
ing v = ±4 trace) and one dip indicate the CR positions. The 
white traces overlaid on the IR data are best Lorentzian fits 
from which we determine the resonance position and width. 
The integrated intensity of the resonances increases roughly 
linearly with B. The widths of the v = ±4 transitions are 
generally twice the widths of the transitions at higher filling 
factors. A typical halfwidth of 14 meV for the lowest LL tran- 
sitions corresponds to a carrier scattering rate of 2 x 10 13 s _1 
and a mobility of 3000 cm 2 /Vs, in reasonable agreement with 
values from the transport data. 

Our complete set of CR energies is shown in Fig. 2, com- 
prising eight separate transitions at five B-fields. The tran- 
sitions are identified by filling factor, and plotted as a func- 
tion of B. The data can be clearly distinguished by their B- 
dependence. While the data in panels (b-d) follow a VB be- 
havior (dashed lines), those in panel (a) are roughly linear in 
B. Additionally, in all panels there is an apparent asymmetry 
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FIG. 1: Representative traces showing normalized IR transmission 
data. In each, the resonance dip is due to the CR transition at filling 
factor v = -8 for B = 18 T, v = 12 for 14 T, and v = 16 for 
10 T, while the resonance peak corresponds to the CR at v = ±4 of 
the normalizing scans. White overlay traces are best Lorentzian fits, 
with arrows indicating dip and peak positions. The traces are offset 
vertically. Upper left inset: quantum Hall effect of bilayer graphene 
sample measured in situ at B = 18 T. Upper right inset: schematic 
of the zero-field dispersion relation of bilayer graphene. 



between electron and hole data, with the energies of equiva- 
lent transitions differing by up to 20%. Such asymmetry, al- 
though on a smaller scale, has also been observed in IR data 
of single layer graphene lflo[[l6ll . Its origin is presently unre- 
solved, but is speculated to be either intrinsic ifioll or caused 
by residual charged impurities (17, 181. 

At first glance the data in Fig. 2 appear to match theoretical 
expectations: for the lowest transitions the behavior is close 
to linear in B, while transitions at higher energies follow VB. 
In particular, we find a good fit of Eq. 2 to both v = ±4 data 
sets in panel (a), where V g set to v = —4 (+4) corresponds 
to an intraband transition n = — 1 — > (0 — > 1). Assuming 
71 = 0.35 eV IU9I1 and leaving c as a fitting parameter we 
deduce c = 1.15 x 10 6 m/s and 1.02 x 10 6 m/s for electrons 
and holes, respectively, shown as solid lines in Fig. 2(a). In 
contrast, for Fig. 2(b-d) reasonable fits can only be achieved 



using a value 71 -c 0.35 eV. The dashed lines in panels (b-d) 
are fits with 71 = 0, leading to E oc VB, which clearly is 
an improvement over using 71 = 0.35 eV. This is surprising, 
since although the spectrum given by Eq. 2 does approach a 
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FIG. 2: LL transition energies as a function of B. Open red (filled blue) symbols represent transitions between electron (hole) LLs. Solid lines 
represent best fits using Eq. 2 with 71 = 0.35 eV and the band velocity c as a fitting parameter. Dashed lines are fits to Eq. 2 using 71 = 0. 



\/B dependence for energies of the order of 71, we see this 
behavior at energies as low as E sa 71 /4. Moreover, one set 
of 71 and c values in Eq. 2 should describe all LL energies. 
Yet, for panels (b-d) it is clear that setting 71 = and fitting c 
over a range of values gives better fits to the data. These fits— 
which mimic the dispersion of single layer graphene- result 
in c = 0.6 — 0.7 x 10 6 m/s for the v = ±8 data, rising to 
0.8 - 0.9 x 10 6 m/s and 0.9 - 1.0 x 10 6 m/s for v = ±12 and 
±16. For clear comparison, the solid lines in Fig. 2 (b-d) show 
fits retaining 71 = 0.35 eV, and give c = 1.0 — 1.2 x 10 6 m/s. 
In panel (a), the dashed line is a fit using 71 = 0, which for 
these two transitions is clearly a poor fit. Therefore, instead 
of showing a smooth transition from a linear-in-B to a linear- 
in-\/B dependence with increasing energy, these data reflect 
a band structure that suddenly becomes linear at energies far 
lower than expected. 

While Fig. 2 attempts to fit the parameters of a theoretical 
dispersion to match our results, our CR data can be presented 
in a fashion that directly arrives at the DOS. The procedure is 
as follows (see inset to Fig. 3): the total number of states in 
a LL is till = deB/h, where d is the LL degeneracy. This 
can be viewed as a condensation of the zero-field DOS in the 
vicinity of each LL energy into an equal number of states in 
the LL; half of them deriving from energies above, the other 
half from energies below the LL. This statement is exact for 
single layer graphene, while for bilayers it is exact in the high 
LL limit, and has no larger than 7% error for LL n = 1. The 
zero-field DOS is then derived from 

,w v d'eB 1 

g(E„,„+i) = -j-Tp , O) 

where AE n! „ + i is the measured transition energy, n and n+ 1 
are the initial and final LL, and n+1 is the energy midway 
between these LLs. The term d'eB/h is the number of states 



in the zero-field DOS that lie between LLs. In bilayers, d' ss 4 
for all transitions except for those to/from n = 0, where we 
approximate it as d' = 6 accounting for the 8 -fold degeneracy 
at Eo = and the small error noted above. E' n n+1 is simply 
determined by "stacking" the transition energies, e.g. E' 23 = 
AE ,i + AEi l2 + 5AE 2 , 3 - 

This method works very well for single layer graphene, us- 
ing previous data for the n = — > 1, — 1 — > 0, 2 — > 3 and 
—3^—2 LL transitions from Ref. [9[], as shown in the lower 
trace of Fig. 3. The so-derived DOS is linear in energy and 
has an intercept near the origin, reproducing closely the form 
of the single particle DOS, g(E) = 2E/(Trh 2 c 2 ). Fitting the 
data with this DOS gives a band velocity c = 1.04 x 10 6 m/s. 

The equivalent data reduction performed on our bilayer data 
appears in the upper left portion of Fig. 3. Although the scat- 
ter of the data is considerably higher than for the single layer 
data, overall, the resulting energy dependence is linear with a 
constant offset, as expected from theory for a hyperbolic dis- 
persion. This is quite satisfactory, since no particular disper- 
sion was assumed to generate the graphs of Fig. 3. As in Fig. 
2 we observe an electron-hole asymmetry. From least square 
fits to the expected DOS, g(E) = (2E + 7i)/(7r^ 2 c 2 ), we es- 
timate 71 = 0.43 eV and 0.52 eV for electrons and holes, 
respectively. These values are close to those previously re- 
ported on bilayer graphene in ARPES and Raman experiments 
|jl2[[l9ll . For both electrons and holes, the band velocity is es- 
timated at c = 1.2 x 10 6 m/s. This is higher than for single 
layer graphene, but includes the single layer value within er- 
ror. In fact, closer inspection of the bilayer data of Fig. 3 
reveals that the scatter is not random, but correlated. 

Two trends stand out in Fig. 3. Since the symbols used here 
are identical to those in Fig. 2, it is apparent that data points 
grouped by transition (y) form lines that are much steeper 
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FIG. 3: DOS as a function of energy, derived from LL transition en- 
ergies for single layer and bilayer graphene. Open red (filled blue) 
symbols represent electron (hole) data. For single layer data OJ, in- 
verted triangles are n = — 1 — > and — > 1 LL transitions; thin 
diamonds are —3 — > —2 and 2 — > 3. For bilayer data, the symbols 
are the same as Fig. 2. The solid line fit through the single layer data 
yieldsc = 1.04 x 10 6 m/s. Dashed (dotted) lines are fits to the bilayer 
electrons (holes) yielding 71 = 0.43 eV and 0.52 eV, respectively. 
For both, c = 1.20 x 10 6 m/s. Symbols overlaid with + ( x ) indicate 
data at B = 10 (16) T. 

than the average slope and tend to a y-intercept (71 /c 2 ) much 
closer to the origin. The increased slopes imply generally 
lower band velocities (by ~ 30%) than the average. These 
trends are consistent with the 71 = fits in Fig. 2 (b-d), in 
which 71 and c were derived from data grouped by transition. 
On the other hand, data grouped by B field (e.g. B= 10 T and 
16 T are marked + and x in Fig. 3) show roughly the same 
slope as the average (1.2 x 10 6 m/s), but are offset in the ver- 
tical direction and extrapolate to different y-intercepts. This 
implies a band velocity not far from the average, and a 71 that 
generally falls between 0.3 eV to 0.7 eV, depending on B. The 
origins of these correlations in Fig. 3 are unclear. However 
they are foreshadowed by our inability to fit the data of Fig. 
2 with one 71 -c parameter set, and are therefore not the result 
of approximations made in arriving at Fig. 3. In any case, in 
spite of these approximations, the DOS in Fig. 3 derived for 
single layer graphene is well behaved. An obvious culprit is 
the naive interpretation of LL transition energies in terms of 
a one particle band structure. LL transitions in single layer 
graphene have already shown inconsistencies with a one par- 
ticle picture at the 5% level, with the deviation attributed to 
many particle effects as large as 30% J9],|2(J. Similar effects 
in bilayer graphene have been antici pate d, as Kohn's theorem 
does not hold for bilayers either Il5ll2lll . 

In summary, we have performed IR absorption measure- 



ments in bilayer graphene and observed multiple cyclotron 
resonance features. We find the B-field dependence of these 
LL transitions to be roughly linear for the lowest transition for 
both electrons and holes, while the remaining higher LL tran- 
sitions display a clear \/B-dependence. This highly unusual, 
varying B-dependence reflects qualitatively the expected be- 
havior for the theoretically predicted hyperbolic dispersion re- 
lation. In detail, however, the transitions do not follow such a 
simple picture and suggest that additional physics contributes 
to the transition energies. 
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